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ABSTRACT 

The radio emitting X-ray binary GRS 1915+105 shows a wide variety of X-ray and 
radio states. We present a decade of monitoring observations, with the RXTE-ASM 
and the Ryle Telescope, in conjunction with high-resolution radio observations using 
MERLIN and the VLBA. Linear polarisation at 1.4 and 1.6 GHz has been spatially 
resolved in the radio jets, on a scale of ~ 150 mas and at flux densities of a few mJy 
Dcpolarisation of the core occurs during radio flaring, associated with the ejection of 
relativistic knots of emission. We have identified the ejection at four epochs of X-ray 
flaring. Assuming no deceleration, proper motions of 16.5 to 27 mas per day have been 
observed, supporting the hypothesis of a varying angle to the line-of-sight per ejection, 
perhaps in a prcccssing jet. 

Key words: ISM: jets and outflows - X-ray binaries: individual (GRS 1915+105). 



1 INTRODUCTION 

GRS 1915+105 has proved to be one of the best labora- 
tories for the study of relativistic physical environments, 
due to its high-brightness, periodic outbursts of superlu- 
minal ejecta and relative proximity. Since the first detec- 
tion in August 1992 with the W ATCH instrument on boar d 
the X-ray telescope GRAN AT (|Castro-Tirado et al.|[l992T ). 
GRS 1915+105 has demonstrated some of the most spec- 
tacular high-energy physics within the Galaxy. Shortly af- 
ter its discov e ry, a radio counterpart was discovered by 
iMirabel et alj (Il993a) and later an in frared component was 
also identified I Mirabel et al.|[l993bl) . Its radio light-curve 
shows a highly variable and complex structure, with spa- 
tially resolved features that can be followed over many days. 

A clear correlation between the X-ray, in- 
frar ed and radio emi ss ion was quickly established 
by IMirabel fc Rodriguez! (|l994f ). with rapid time vari- 
ability in each band. Radio flaring was found to correspond 
to a rapid change in the hard X-rays and possible production 
of high-energy gamma rays. In March 1994, 20 cm VLA ob- 
servations of GRS 1915+105 d etected the first superlumina l 
motion of a Galactic source (|Mirabel fc Rodriguez! 1 1994(1 . 
This major breakthrough provided the direct evidence 
of relativistic jets and an extreme physical environment 
within the Galaxy. The name "microquasar" was coined 
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due to their obvious similarities with their extra-galactic 
counterparts, quasars. 

The launch of the Rossi X-ray Timing Explorer (RXTE) 
satellite, in December 1995, signified the start of a long-term 
monitoring campaign of X-ray binaries. The All Sky Moni- 
tor (ASM) on board the RXTE ha s taken daily observat ions 
of GRS 1915+105 since its launch. iGreiner et alj (|l996l ) de- 
tected unusual X-ray variability on time scales of under one 
second to days. The RXTE-ASM lightcurve was found to be 
both highly complex and structured, due to instabilities in 
the accretion disc. Detailed observations with t he RXTE's 
Proportional Counter Ar ray (PCA) instrument (|Chen et al.l 
ll997tlBelloni et al.lll997l ) identified different spectral states, 
including a low-hard state dominated by a power-law and a 
high-soft state with a strong disc-blackbody component. A 
transition between such states is believed to be associated 
with the ejection of superluminal plasmons (jFender et al.l 
1 1999(1 . 

Radio polarisati on was first detecte d with MERLIN 
in GRS 1915+105 by iFender et all (|l999l ) in the first four 
epochs (i.e. first four days) of the October 1997 flare at a 
frequency of 5 GHz. A clear asymmetry in the location of 
the polarised emission was found, showing the strongest de- 
tection of linearly polarised emission within the approaching 
jet, weaker in the receding jet and no detection (< 3% of the 
total flux density) within the core. The polarisation position 
angle (PA) was found to rotate by at least 75° over the four 
days, leading to the suggestion that the changes in polarisa- 
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Epoch (MJD) 


Frequency / Array 
(GHz) 


(mjy) 


LP total 
(mJy) 


Fraction 

(%) 


Noise 
(/ijy/beam) 


Spectral* 
Index 


Resolved? 


*1997 October 31 (50752) -core 


2.27 / VLBA 


105 ± 2 


< o.i* 


< 0.5* 


430 






*1997 October 31 (50752) -jet 


2.27 / VLBA 


226 ± 2 


5.6 ± 0.3 


2.5 


430 




YES 


*1997 October 31 (50752) -core 


8.3 / VLBA 


33 ±1 


< o.it 


< 0.5* 


460 


- 


*1997 October 31 (50752) -jet 


8.3 / VLBA 


65 ± 3 


4.7 ± 0.3 


7.2 


460 






1998 June 9 (50973) 


1.6 / MERLIN 


145.1 ± 0.4 


9.5 ± 0.3 


6.9 ±0.2 


300 


-0.3 


NO 


1999 April 2 (51270) 


1.6 / MERLIN 


238.4±0.3 


3.5 ± 0.3 


1.5 ±0.2 


530 


-0.6 


NO 


1999 November 18 (51500) 


1.6 / MERLIN 


221.3±0.8 


6.9 ±0.3 


3.2 ±0.2 


860 


-0.4 


NO 


1999 November 22 (51504) 


1.6 / MERLIN 


37.1±0.7 


sC 0.1+ 


< 0.4* 


300 


-0.2 


NO 


1999 December 28 (51540) 


1.6 / MERLIN 


89.1±0.7 


12.7 ± 1.6 


14.1 ±0.2 


400 


-0.4 


NO 


2003 March 6 (52704) 


1.6 / MERLIN 


31.3 ±0.3 


1.6 ±0.3 


5.2 ± 1.0 


240 


-0.4 


NO 


2003 March 25 (52723) 


1.6 / MERLIN 


125.8 ±0.1 


7.1 ± 0.2 


5.6 ±0.2 


140 


0.2 


NO 


2003 April 18 (52747) -core 


1.6 / MERLIN 


135.4 ±0.4 


3.5 ± 0.2 


2.6 ±0.2 


270 


0.0 


YES 


2003 April 18 (52747) -jet 


1.6 / MERLIN 


20.0 ±0.3 


1.3 ±0.1 


6.5 ±0.5 


270 


2003 May 09 (52768) 


1.6 / MERLIN 


43.0 ±0.1 


2.7±0.1 


6.3 ±0.2 


160 


0.2 


NO 


2003 June 15 (52805) -core 
2003 June 15 (52805) -jet 


1.6 / MERLIN 
1.6 / MERLIN 


41.4 ±0.1 

67.5 ±0.1 


SC 0.1* 
7.6 ±0.4 


sC 0.2* 
11.9 ±0.6 


390 
390 


-0.4 


YES 


2006 December 24 (54093) 
2006 December 24 (54093) 


1.4 / MERLIN 
1.6 / MERLIN 


145.0 ±0.5 
156.9 ± 1.2 


7.7 ± 0.6 
11.5 ±0.5 


5.3 ±0.4 
7.3 ±0.5 


380 
380 


-0.1 


NO 


2006 December 27 (54096) 
2006 December 27 (54096) 


1.4 / MERLIN 
1.6 / MERLIN 


20.7 ±0.5 
27.0 ±0.5 


^0.1* 
SC 0.1 1 " 


^ 0.5* 
«: 0.5* 


230 
230 


-0.4 


NO 


2006 December 28 (54097) 


1.4 / MERLIN 


25.8 ±0.4 


< 0.1 T 


^ 0.4* 


210 




NO 


2006 December 28 (54097) 


1.6 / MERLIN 


26.2 ±0.4 


^ o.it 


^ 0.4* 


210 




2007 January 04 (54104) 
2007 January 04 (54104) 


1.4 / MERLIN 
1.6 / MERLIN 


21.0 ±0.5 
22.2 ±0.5 


<0.1* 

s; o.it 


0.5* 
< 0.5* 


260 
260 


-0.3 


NO 



Table 1. Journal of observations and results: Date of observation, observed frequency and array, total flux, total linearly polarised 
flux, fraction of total flux that is lin early polarised, syste matic noise, spectral index and detection of resolved structure. [Notes: * The 
observations were first published by iDhawan et al.ll2000l: * Estimates are based on a upper limit of detection; * The spectral index is 
derived from the total MERLIN flux (i.e. core + jet) and the 15 GHz RT data, where S oc v a .] 



tion could be due to Faraday effects (i.e. changes in the Fara- 
day depth as the plasmons move away from the core) , imply- 
ing a change in rotation measure of ARM > 300 rad m -2 . 

This paper describes observations taken over a decade 
with MERLIN at 18 and 21 cm. The observations were ei- 
ther in conjunction with INT EGRAL observa t ions as part 
of th e Galactic Plane Survey (|Bird et al J 12004 iFuchs et all 
|2004 >. or were triggered by flare events, found by other radio 
telescopes. Polarisation behaviour, structural variations and 
the relationship of activity in the radio regime to the X-ray 
behaviour were investigated. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 RXTE All Sky Monitor 

The All Sky Monitor (ASM) instrument on-board the RXTE 
has been monitoring the sky since March 1996 and the data 
presented here covers the period from March 1997 to 2007. 
With each orbit of the RXTE, the ASM surveyed ~ 80% 
of the sky to a depth of 20 — 100 mCrab, making approxi- 
mately ten observations of a source per day. A more detailed 
descri ption of the RXTE- ASM can be found in lLevine et alj 
(1 19961 ). 

Each individual pointing, or dwell, was a 90 second in- 
tegration of the source, with intensities measured in three 
energy bands of 1.5 — 3, 3 — 5, and 5 — 12 keV. The Crab 
Nebula flux between 1.5 — 12 keV corresponds to about 75 
ASM counts s _1 . To calculate the spectral hardness (HR2), 
individual dwells were averaged into daily points and the 



ratios between the 5 — 12 keV and 1.5 — 3 keV energy bands 
were taken. 



2.2 Ryle Telescope 15 GHz monitoring 

The Ryle Telescope (RT) is a linear east-west radio interfer- 
ometer located at the Mullard Radio Astronomy Observa- 
tory, UK. The array operates at a frequency of 15 GHz with 
associated baselines between approximately 18 metres and 
4.8 kilometres. 

An extensive monitoring campaign began with the RT 
in 1996 of a few radio-bright X-ray binaries (including 
GRS 1915+105), which coincided with the launch of the 
RXTE satellite. Observations of target sources were inter- 
leaved with a nearby phase calibrator (B1920+154) and the 
flux-density scale was set by short scans of 3C 48 and 3C 286. 
The data were sampled every eight seconds and averaged 
into five minute data points with an RMS of ~ 2 mjy. 

iPoolev fc Fenderl l] 19971) describe the details of this pro- 
gramme, indicating the detection of a 20 - 40 minute quasi- 
periodic variation of GRS 1915±105 at a frequency of 15 
GHz, possibly coupled with variations in the soft X-rays. 
The results presented in Section [3] were taken during moni- 
toring periods with the RT and RXTE- ASM. Figures [1] and 
show data collected during 1997, 1998, 1999, 2003 and 
2006 - 2007, each for a period of 365 days (note RT data 
was not available in 2006 - 2007). 
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Figure 1. X-ray observations with the All Sky Monitoring instrument on board the RXTE satellite and 15 GHz radio monitoring with 
the Ryle Telescope in 1997 and 1998. The dotted line denotes a MERLIN 18 cm epoch (1998 June 9), showing th at the source was in th e 
'plateau-state'. [N.B. The solid line in 1997 marks a VLBA observation (1997 October 31) originally published bv lDhawan et all J2000T> ] ■ 



2.3 RATAN-600 

The RATAN-600 radio telescope consists of a 576 metre cir- 
cle of radio reflectors, located at the Special Astrophysical 
Observatory, Russia. The telescope can operate at multiple 
radio frequencies between 610 MHz and 30 GHz. 

RATAN-600 observations were carried out as part 
of a lon g-term monitori n g programme on var iable X-ray 
sources (jTrushkinl |2000| ; frrushkin et"aH |2006( ). Observa- 
tions of GRS1915+105 were taken at TO, 2.3, 4.8, 7.7 and 
11.2 GHz during the 2006 - 2007 flare: on 2006 December 
22 (MJD 54091), 24 (MJD 54093), 27 (MJD 54096) and 



2007 January 04 (MJD 54104). A standard continuum radio 
system with cryogenic receivers made observations at 4.8, 
7.7 and 11.2 GHz. At lower frequencies, a low-noise HEMT- 
based radiometer operated at 1 and 2.3 GHz. The observa- 
tions were carried out with the 'Northern sector' antenna of 
RATAN-60(E 



http: / /w0. sao.ru/ratan 
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Figure 2. X-ray observations with the All Sky Monitoring instrument on board the RXTE satellite and 15 GHz radio monitoring with 
the Ryle Telescope in 1999, 2003 and 2006 - 2007 (no RT data were available in 2006 - 2007). The dotted lines denote the MERLIN 
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2.4 VLBA 

The Very Long Baseline Array (VLBA) is a network of 
remotely controlled radio telescopes distributed across the 
USA. The Maximum baselines are over 8500 kilometres pro- 
ducing a sub-milliarcsecond resolution at centimetre wave- 
lengths. 

Observations taken on 1997 October 31 (MJD 50752) 
with the VLBA have been presented by iDhawan et al.l 
120001 ). The VLBA data have been reprocessed in order to 
study linearly polarised emission from the source. The ob- 
servations were originally triggered by daily monitoring with 
the Green Bank Interferometer (GBI) during a period of ra- 
dio flaring. It was possible to obtain data from multiple fre- 
quencies by means of a frequency-selective optics allowing 
the simultaneous use of 3.6 and 13 cm (8.4 and 2.27 GHz) 
receivers. 



2.5 MERLIN 

The Multi-Element Radio Linked Interferometer Network 
(MERLIN) array connects a total of seven telescope across 
the UK: Mark II and Lovell at Jodrell Bank, the 32 me- 
tre at Cambridge and 25 metre antennas at Knocking, 
Darnhall, Pickmere and Defford (although the Lovell was 
not included in all observations). MERLIN observations of 
GRS 1915+105 have been made during 14 epochs since 1998 
(see Table [TJ at 1.4 and 1.6 GHz, where the angular resolu- 
tion was 150 and 130 mas, respectively. Note t hat MERLIN 
observ atio ns at 5 GHz have been des cribed bv lFender et al.l 
|l99i) and iMiller- Jones et all (|2007n . 

The antennas have orthogonal right (R) and left (L) 
circular polarisation feeds from which RR, LL, RL and LR 
cross-correlations were formed, measuring 15 MHz band- 
widths correlated into 1 MHz channels. Data correlation was 
performed by the dedicated MERLIN correlator located at 
Jodrell Bank, UK. Preliminary flux calibration and gain ele- 
vation corrections were applied to the data using local MER- 
LIN software and transformed into fits format. The data 
were loaded into NRAO's Astronomical Image Processing 
System (aips) and bad data were removed using the task 
ibled. The flux density scale was determined from observa- 
tions of 3C 286, using the three short baselines between the 
Mark II, Pickmere and Darnhall telescopes. A flux density 
of 13 .5 Jy at 1.6 GHz was assumed for 3C 286 (|Baars et al.l 
1977), from which we derived a flux density for the point 
source calibrator OQ 208. Scans of GRS 1915+105 were in- 
terleaved with the phase reference source, 1919+086, with 
a 5:3 minute cycle time respectively, allowing the use of an- 
tenna amplitude and phase-gain solutions from the calibra- 
tor source using the tasks CALIB and CLCAL. 

The AIPS task PCAL was used on 1919+086 to determine 
the effective feed polarisation parameters for each individual 
telescope and these were written into the antenna file (or AN 
table) for future correction. This corrects the leakage of flux 
from one cross-polarised mode to the other, assuming the 
phase and amplitude of the source has been correctly cal- 
culated and is flat across the band (i.e. phase referenced). 
The instrumental response due to changes in the parallac- 
tic angle was also removed using observations of the phase 
reference source 1919+086. 

Next, any systematic phase delay between the right and 
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Figure 3. RATAN observations of the 2006 - 2007 flare at the 
frequencies 1.0, 2.3, 4.8, 7.7 and 11.2 GHz. 



left hands needed to be corrected by a source with a known 
polarisation angle. A phase delay difference between right 
and left hands was removed before performing PCAL. Whilst 
applying the solutions from PCAL, a rotation of Q + iU 
corrected the R-L phase delay was achieved using the AIPS 
task CLCOR. The known position angle of the polarization of 
3C 286 was used to calibrate the phase differences between 
L and R for each antenna. 

Finally the target source went through a few rounds 
of phase-only self-calibration. The calibrated uv data of 
GRS 1915+105 were then Fourier transformed and the 
CLEAN algorithm was applied using the AIPS task IMAGR. 
Each image was made using approximately 12 hours of 
MERLIN data. Table [T] shows the total flux density and the 
polarized flux density from the images. Also shown are the 
results of VLBA observat i ons w hich were re-processed using 
data from lDhawan et alj (|2000f ). These values are presented 
for a multi-wavelength comparison of the polarisation levels 
in the core and jet. 



3 RESULTS AND POLARISATION 
MEASUREMENTS 

The target source was detected in all 14 MERLIN obser- 
vations, with a total flux density varying between 20 and 
300 mjy at 1.4 or 1.6 GHz. Most epochs revealed a to- 
tally unresolved core, except six which showed a slightly 
extended structure with a position angle to the south-east. 
Linear polarisation was found with MERLIN in 10 out of the 
14 epochs and at both frequencies of the VLBA observation. 
Individual 1 MHz MERLIN channels were able to detect the 
linear polarisation with a position angle accuracy of ±10°; 
however, no significant change in the position angle of the 
polarisation was found across either bandpass (i.e. between 
1401 - 1415 MHz or 1651 - 1665 MHz). Likewise a clear 
change in the position angle between 1.4 and 1.6 GHz was 
inconclusive due to high variability and low image fidelity. 

Continual monitoring with the RT has shown frequent, 
variable radio flaring over the last decade. Strong radio flares 
(> 100 mjy) are associated with a hardening of the X-ray 
emission (HR2 > 1.4) and a persistently soft flux (2—12 keV) 
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Figure 4. MERLIN images of GRS 1915+105 at 1.6 GHz (top left to bottom right) centred on Right Ascension 19 15 11.54801 
Declination +10 56 44.7662, on 1998 June 3, 1999 April 02, November 18, 22, December 28 and 2003 March 06. The restoring beam is 
shown in the lower corner of each image and the stick vectors represent the mean position angle of linear polarisation with an intensity 
of 12.5 mjy per arcsecond. The contours are 3cr X — 1, 1,2, 4, 8, 16....n. (a) - 1998 June 09. 3cr = 0.90 mjy with a peak flux density = 
95.1 mjy. (b) - 1999 April 02. 3ct = 1.6 mjy with a peak flux density = 221.7 mjy. (c) - 1999 November 18. 3ct = 2.6 mjy with a peak 
flux density = 194.8 mjy. (d) - 1999 November 22. 3<r = 0.9 mjy with a peak flux density = 8.0 mjy. (e) - 1999 December 28. 3tr = 1.2 
mjy with a peak flux density = 49.2 mjy. (f) - 2003 March 06. 3a" = 0.733 mjy with a peak flux density = 22.3 mjy. 



High-resolution radio observations of GRS 1915+105 7 




ARC SEC ARC SEC 
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Figure 5. MERLIN images of GRS 1915+105 at 1.6 GHz (top left to bottom right) centred on Right Ascension 19 15 11.54801 
Declination +10 56 44.7662, on 2003 March 24, April 18, May 09, June 15, 2006 December 24 and 27. The restoring beam is shown in 
the lower corner of each image and the stick vectors represent the mean position angle of linear polarisation with an intensity of 12.5 
mjy per arcsecond. The contours are 'Aa X — 1, 1, 2, 4, 8, 16. ...n. (a) - 2003 March 25. 3a = 0.422 mjy with a peak flux density = 110 mjy. 
(b) - 2003 April 18. 3a = 0.822 mjy with a peak flux density = 86.5 mjy. (c) - 2003 May 09. 3a = 0.48 mjy with a peak flux density = 
43.0 mjy. (d) - 2003 June 15. 3a = 1.18 mjy with a peak flux density = 49.7 mjy. (e) - 2006 December 24. 3a = 1.14 with a peak flux 
density = 144.5 mjy. (f) - 2006 December 27. 3a = 0.69 mjy with a peak flux density = 17.7 mjy. 
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Figure 6. MERLIN images of GRS 1915+105 at 1.6 GHz centred on Right Ascension 19 15 11.54801 Declination +10 56 44.7662,on 
2006 December 28 and 2007 January 04. The restoring beam is shown in the lower corner of each image and the stick vectors represent 
the mean position angle of linear polarisation with an intensity of 12.5 mjy per arcsecond. The contours are 3u X -1, 1, 2, 4, 8, 16. ...n. 
(a) - 2006 December 28. 3cr = 0.64 mjy with a peak flux density = 19.8 mjy. (b) - 2007 January 04. 3<r = 0.78 mjy with a peak flux 
density = 14.9 mjy. 



of fn 50 counts/second with a low RMS. This state has pre- 
viousl y been identified as the 'plateau' state bv lFender et al.l 
|l999l ). Weaker radio flares (^C 100 mjy) do not clearly cou- 
ple to the X-ray emission, suggesting an underlying mecha- 
nism not related to X-ray state changes. 



Individual events 

October 1997 (MJD 50730-50750) 

In October 1997 the RT observed the first period of a sus- 
tained radio outburst, begining around 1997 October 04 
(MJD 50725) and lasting approximately 30 days, ending 
with a short flare peaking at 220 mjy (at 15 GHz). A 
series of 5 GHz M ERLIN observations were reported by 
l|Fender et al.|[l999h , showing t he first Galacti c supe rluminal 
proper motion with MERLIN. iDhawan et al.l (|2000h also re- 
ported near simultaneous VLB A observations, showing an 
extended jet to the south-east (SE). Our analysis of the 
VLBA data has detected linearly polarised emission from 
within the extended jet at 2.5% and 7.2% fraction of the to- 
tal intensity (Lotai) at frequencies of 2.27 GHz and 8.3 GHz, 
respectively, but none within the core. 

Over this period the RXTE-ASM detected a sudden 
quenching of the soft X-rays, with a sharp rise in the hard- 
ness ratio, HR2 > 1.6. The ASM count rate then remained at 
a level of ss 50 counts/second, with a low RMS fluctuation. 
This state ended with a transition to a different state with 
a short radio and soft X-ray flare of w 150 counts/second, 
as sociated with t he eje ction of a superluminal knot reported 
in lFender et~ai1 (|l999h . 

April - July 1998 (MJD 50905-51020) 

A major radio outburst was detected by the RT start- 
ing around 1998 April 04 (MJD 50907), which lasted 



for « 110 days. During this period the RT also detected mul- 
tiple short flares, with the brightest flare reaching 362 mjy 
(at 15 GHz). The X-ray emission returned to the plateau 
state of w 50 counts /second, as seen in October 1997, and 
was notably harder during this outburst, with the ASM 
hardness ratio reaching HR2 ~ 2.0. The multiple short radio 
flares can also be seen to be associated with the soft X-ray 
flares throughout this quenched state (as shown in Figure[l|. 

A MERLIN observation on 1998 June 09 (MJD 50973) 
was made six days after one of the short flares (Figure |4]-a) , 
showing an extended jet to the SE of the core. A total flux 
of 145 mjy (at 1.6 GHz) was detected compared to a RT 
observation of 79 mjy (at 15 GHz), giving a spectral index of 
a — —0.3 (where S oc u a ). Linear polarisation was detected 
by MERLIN at a total flux of 9.5 mjy, which was ~ 7% 

Of /total- 



April 1999 (MJD 51270) 

Paradoxically, the strongest radio flare detected in these 
MERLIN observations occured during a period of relative in- 
activity in radio and X-rays. On 1999 April 02 (MJD 51270) 
MERLIN 1.6 GHz observations detected a peak flux density 
of nearly « 300 mjy, which quickly decayed to m 170 mjy. 
Little evidence of a jet was found in the image (Figure 2}b) 
with only a marginally resolved core detected. Linearly po- 
larised emission was detected in the core at a total flux of 
3.5 mjy, which was ~ 1.5% of /total- 

The RT did not detect a sustained period of radio out- 
burst and only measured a short flare with a relatively mod- 
est flux density of ~ 60 mjy (at 15 GHz), giving an approx- 
imate spectral index of a ~ —0.6. RXTE-ASM observations 
during April 1999 also had a relatively low count rate of 
45 — 60 counts/second, showing no sign of an X-ray flare 
during this period. 
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November - December 1999 (MJD 51500-51550) 

The latter two months of 1999 had a modest period of 
radio activity following a period of relative quiet. The 
RT measured short flares which peaked at 126 mjy and 
171 mjy (at 15 GHz), and a sustained period of radio out- 
bursting occurred for ~ 15 days. Although continual RXTE- 
ASM observations were not available, X-ray counts peaked 
at « 210 counts/second and the hardness radio reached 
HR2 ~ 2. 

MERLIN observations on 1999 November 18 
(MJD 51500) detected another short, strong radio flare of 
221 mjy (at 1.6 GHz), whilst the RT measured ~ 90 mjy 
(at 15 GHz), making a spectral index of a = —0.4. Again 
only a marginally resolved core was found (Figure |4]-c), with 
linearly polarised emission at a fractional level of ~ 3.2%. 
Four days later on 1999 November 22 (MJD 51504) MER- 
LIN observed that the radio flare had quickly dissipated to 
37 mjy (Figure |Ud), showing no linearly polarised emission 
(^ 0.4%). At 15 GHz there was a similar flux density of 
30 mjy, resulting in dramatic change in the spectral index, 
which reduced to a = —0.2. 

Finally on 1999 December 28 (MJD 51540) a small flare 
of 90 mjy was detected with MERLIN at 1.6 GHz and RT 
at 15 GHz with 40 mjy (a ~ —0.4). A very unusual struc- 
ture was seen to north of the core (FigureUJe), not normally 
associated with the position angle of previous jets. The lin- 
ear polarisation was also found to be unusually high for a 
frequency of 1.6 GHz with polarised flux of 12.7 mjy, i.e. 

« 14% Of /total- 



February - July 2003 (MJD 52695-52825) 

The most spectacular event occured in the first half of 2003, 
in a giant outburst of radio activity, with multiple short 
radio and X-ray flares. The event began on 2003 March 02 
(MJD 52700) with the RT observing a short flare at 15 GHz 
of 151 mjy, followed 14 days later with peak flare of 372 mjy 
and then a slow decay over the next ~ 100 days. The RXTE- 
ASM, as seen in 1997, observed a quenching of the soft X- 
rays, with a strong rise in the hardness ratio and the emission 
remaining at a level of ~ 50 counts/second, with short X-ray 
flares occuring with radio flares. 

Five MERLIN observations were taken over this period 
of activity [starting four days after the initial radio flare on 
2003 March 06 (MJD 52704)]. This epoch showed only a 
slight NW-SE extension (Figure |3rf), suggesting a small jet 
associated with the initial outburst with a weak flux density 
of ~ 31 mjy (at 1.6 GHz), whilst the RT measured 15 mjy 
(at 15 GHz), giving a spectral index of a = —0.4. Only weak 
linearly polarised emission was found near the small jet, at 
a fraction of the total emission of 5.2%. 

The second observation occured on 2003 March 25 
(MJD 52723), seven days after the peak flare was detected 
with the RT. The MERLIN map showed very little exten- 
sion (Figure [SJ-a), but had a variable core of 120 mjy (RT 
data was not available). Polarised emission was again de- 
tected within the core at a fractional level of ~ 5.6 %. No 
X-ray flares were detected with the ASM instrument over 
this time. 

On 2003 March 05 (MJD 52734), 35 days after the ini- 
tial radio flare, the RXTE detected a short, but relatively 



strong and soft X-ray flare of 150 counts/second. This was 
followed by another increase in radio emission and then a re- 
turn to the constant hard X-ray state. MERLIN (1.6 GHz) 
observations 13 days after this short event, on 2003 April 18 
(MJD 52747), revealed an extended component of 20 mjy 
to the south-east, with the core at a peak flux density 
of 135 mjy (Figure [5]-b). Measurements at 15 GHz with 
the RT showed a strong flare of 143 mjy, suggesting the 
core had an approximately flat spectrum. Linearly polarised 
emission was detected both within the core (2.4 mjy) and 
the jet (1.3 mjy) at a fractional total emission of 2.6% and 
6.5% respectively. On 2003 May 09 (MJD 52768), 34 days 
after the strong X-ray flare, MERLIN observations showed 
the extended emission had disappeared (Figure [5} c) and an 
unresolved core had decayed to 43 mjy at 1.6 GHz with 
linear polarisation at 6.3% (no RT data were available). 

The 2003 'plateau' state ended on 2003 May 29 
(MJD 52788) with an X-ray flare of 163 counts/second and 
the start of a new period of highly variable soft emission. 
Shortly following this state change, MERLIN observations 
were taken on 2003 June 15 (MJD 52805). The image showed 
an extended jet with a total flux of 67 mjy, which was now 
stronger than the core which had a peak flux 41 mjy (Fig- 
ure (5}d). Furthermore, the linearly polarised emission was 
no longer detectable from within the core and could only be 
detected within the jet (RT was also not available). 



December 2006 - January 2007 (MJD 54091-54105) 

RATAN observations on 2006 December 22 (MJD 54091) 
detected a flare of 350 mjy at 1.0 GHz, 150 mjy at 2 - 8 
GHz and 80 mjy at 11.2 GHz (i.e. a steep spectrum), which 
triggered the final four MERLIN epochs in December 2006 
and January 2007. Observations were taken simultaneously 
at 1.4 and 1.6 GHz by switching between frequencies every 
hour; however, only images produced at 1.6 GHz are shown 
in Figures [5] and [6] as there were no structal differences be- 
tween respective frequencies (i.e. only the peak density flux 
changed). 

The MERLIN observations on 2006 December 24 
(MJD 54093), confirmed the bright flare of ~ 145 mjy and 
~ 157 mjy at 1.4 GHz and 1.6 GHz respectively (Figure [5j 
c) . An unresolved core showed linear polarisation associ- 
ated at a fractional level of 5.3% and 7.3% at 1.4 GHz and 
1.6 GHz respectively. The RATAN spectrum showed that 
the spectrum had flattened, suggesting the optical depth 
had increased. 

However, this was only a short lived flare as on 2006 
December 27 (MJD 54096), just five days after the initial 
detection, the flare had decayed to 21 mjy at 1.4 GHz and 
27 mjy at 1.6 GHz (Figure 0-f) , with RATAN also showing 
the source to be much weaker. A slight extension to the SE 
was detected, suggesting the formation of a weak jet. The 
source remained in a similar weak state for the remaining 
two epochs, on 2006 December 28 (MJD 54097) and 2007 
January 04 (MJD 54104) and all linear emission had reduced 
to < 0.5% (see Figures |6]-a and[6]-b). 

RXTE observations (Figure [2] bottom) showed no obvi- 
ous spectral transition and the radio flare was not associated 
with an X-ray flare. 
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4 DISCUSSION 

This long-term comparative study between the X-ray and 
radio emission from GRS 1915+105 shows a complex and 
varied display of the emission mechanisms in a relatively 
bright X-ray binary. When studied over a ten year period, 
characteristics in both the radio and X-ray have repeated in 
quasi-periodic patterns. 

Broadly summarizing, the observed radio emission is 
attributed to relativistic jets, which are coupled to the ac- 
cretion disk in a poorly understood process. Emission from 
the accretion disk is dominated by X-ray emissions that can 
be modeled by a multi-temperature blackbody and a high- 
energy, non-thermal component. This non-thermal compo- 
nent, fitted with a power-law, is very likely to be due to the 
comptonization of cold electrons surrounding the jet and ac- 
cretion disk. 

Understanding the interplay between the radio jet, 
thermal X-rays and the non-thermal component, requires 
the classification of the different 'states' as described in 
iFoster et all i| 19961 ). Monitoring with the RT and the RXTE- 
ASM, has therefore been used to classify the state of 
GRS 1915+105 into the following: Weak radio emission (<C 
100 mjy) with short-lived flares that do not clearly couple to 
the X-ray emission - this has been identified as the so-called 
'radio-quiet state'; strong radio outbursts (> 100 mjy), typ- 
ically lasting for weeks, that are associated with a hardening 
of the X-ray emission and a persistently soft flux - this has 
been identified as the so-called 'plateau state'; transitions 
from the 'plateau' state accompanied by a short X-ray flare 
and the ejection of superluminal knots - this has been iden- 
tified as the so-called 'flaring state'. 



4.1 No detected jet - 'Radio quiet state/Short 
flares' 

The behaviour of the jet is least understood in the radio 
quiet state of GRS 1915+105. The radio emission is char- 
acterized by a weak flux that is typically less than 5 mjy 
at 15 GHz. Short radio flares are occasionally interleaved, 
which last less than a day and can reach a peak flux of 
100 — 200 mjy at 15 GHz. X-ray emission shows no clear cor- 
relation with the radio during this state and typical RXTE- 
ASM peaks of 100 — 150 counts/second are observed. 

Due to selection effects that limit most observations to 
periods of strong radio emission, only a few high-resolution 
observations have been taken in this state. e -EVN observa- 
tions taken in this state during April 2006 (|Rushton et alj 
12007ft showed n o evidence for any compact jet, unlike the 
jet observed by iFuchs et al.l (|2003T ). suggesting no outflow 
of material. Likewise, MERLIN observations taken in 1999 
showed no large scale structure, despite X-rays flares, sug- 
gesting a de-coupling of the accretion disk to any radio emis- 
sion. 

MERLIN observations taken at the end of 2006, ob- 
served a short radio flare of ~ 150 mjy at 1.6 GHz. This 
emission quickly decayed to less than 20 mjy, again showing 
no large scale ejections. RATAN observations also followed 
the short flare and measured an optically thin spectral in- 
dex that quickly flattened as the flare decayed. The change 
in spectral index indicates that material did flow out of the 
X-ray binary; however, no superluminal knots were formed 



and nor was a steady jet. This is possibly either due to the 
jet not becoming energetic enough or interactions with the 
stellar system preventing the release of material. 



4.2 Steady compact jet - 'Plateau state' 

The plateau state in GRS 1915+105 is analogous to the 
low/hard state of most XRBs, like A0620-00, as the X-ray 
spectrum becomes sli ghtly harder a nd is associated with 
strong radio emission. (Belloni et ail (|200Ch identified these 
X-ray states as class \ using the RXTE-PCA. However, this 
analogy (with the low/hard state) is not entirely accurate as 
the X-ray spectra still rem ains relatively soft compared t o 
other low/hard state XRBs (|McClintock fc Remillardll2006l '). 
Physically this means that the X-ray emission still has a sig- 
nificant contribution from the thermal accretion disk, whilst 
producing steady compact radio jets. As such, this accretion 
state in GRS 1915+105 shows unique coupling between the 
in-falling material and the ejected jets. 

All long periods of radio outburst (typically lasting 
~ 100 days) are asso ciated with the pl ateau state. VLBI 
observation s taken by Fuchs et al.l (120031 ) in April 2003 and 
similarly in iDhawan et al.1 <|2000h showed a steady, compact 
jet on AU-scales. Radio emission was optically thick as the 
jet was partially self-absorbed. This radio flux appeared to 
have a close relationship to the X-rays, which also emitted at 
a steady rate throughout this state. MERLIN observations 
taken in 1998 and over the April 2003 outburst, showed no 
large scale ejections, unless there was an X-ray flare; there- 
fore, this radio emission is constrained relatively close to the 
accretion disk and this is the best state to study the coupling 
between the radio and X-ray emission. 

The cause of the emission in the plateau state is pos- 
sibly due the collapse of an inner accretion disk through 
instabilities in the accretion process. This would cause the 
accretion to be dominated by an advection-dominated pro- 
cess (e.g. ADAF), rather than a standard accretion disk. 
Advection-flow accretion also causes the comptonization of 
a surrounding corona of cold electrons, hence the hardening 
of the X-ray spectra from a more dominant non-thermal X- 
ray component. The inner disk would then be re-filled from 
a surrounding outer disk, and form the compact jet. 



4.3 Ejection of knots - 'Flaring state' 

It has been well established that a state transition in 
GRS 1915+1 05 is associated wi th the ejection of superlu- 
minal knots (|Fender et al.|[l999T ). However, the precise pro- 
genitor of the knot has not been clearly identified nor has 
the exact cause of their formation. iMiller- Jones et al.l (|2005T l 
suggested that an increase in the jet velocity causes internal 
shocks to form and produce the ejected knots. 



Proper motion of knots 

By identifying the precise times of the formation of knots, it 
is then possible to derive their proper mot ions. This assumes 
no dec eleration of the jet, as supported bv lMiller- Jones et al.l 
(|2007l ) who found no deceleration beyond 70 mas from the 
core. The proper motion of an approaching, fj, app , or reced- 
ing, firec, plasmon is given as 
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MERLIN epoch 


Separation 


X-ray flare 


Proper motion 


Intrinsic velocity 


(MJD) 


(mas) 


(MJD) 


(mas d — 1 ) 


(c) 


1998 June 09 (MJD 50973) 


165 ± 15 


50967 


27 ±3 


1.1 ± 0.2 


2003 March 06 (MJD 52704) 


180 ± 20 


52694 


18 ±2 


0.79 ± 0.13 


2003 April 18 (MJD 52747) 


235 ±5 


52733 


16.5 ± 1 


0.74 ± 0.063 


2003 June 15 (MJD 52805) 


305 ±5 


52788 


17.5 ± 1 


0.77 ± 0.063 






52792 


23.5 ±1 


0.97 ±0.07 



Table 2. The proper motions of different ejecta from GRS 1915+105, based on the time after a short X-ray flare. The intrinsic velocity 
is modelled assuming a distance of 12.5 ± 1.5 kpc and flcosd = 0.41. 



/'rt/j 



j3 sin ( 



(1 t/3cos<9) D ' 



(1) 



where 9 is the jet angle to the line of sight, f3 (= v/c) the 
intrinsic jet velocity and D the distance to the source. The 
quantity (3 cos 9 can be measured independently of distance 
from previous observations using 



f3 cos 9 = 



(2) 



where both the approaching and reced ing componen t s were 
observed. MERLIN observations by IFender et ail (|l999l ) 
found f3cos9 = 0.4 1 ± .02 and VLA observations by 
iMirabel fc Rodriguez! (|l994l ) found (3cos9 = 0.323 ±0.011. 

Then, by using Equations [T] and [2] we can find the in- 
trinsic velocity of the jet 



fJ- 2 apP D 2 (l -13 cos 9f 



+ (/3 cos (9) 2 



(3) 



assuming there is no precession of the jet (i.e. (3 cos 9 is a con- 
stant). The MERLIN results presented in Section [3] showed 
four epochs with extended knots that could be identified 
with X-ray flares. Table [2] shows the derived proper mo- 
tions and modelled intrinsic velocities. Assuming a distance 
of 12.5 ± 1.5 kpc from HI absorption measurements, these 
results provide evidence that the velocity of the ejected com- 
ponents was intrinsically different for each event. 

In 1998 June 09 (MJD 50973) a very high proper mo- 
tion was observed of 27 mas d" 1 [c.f. VLA /i apP ~ 17.6 ± 
0.4 mas llMirabel fc Rodriguez! ll994T l; MERLIN n apP ~ 
23.6±0.5 mas ( Fender et al.lll999l )]; this implies an apparent 



motion of 1.5 c (again assuming a distance of 12.5 ±1.5 kpc) 
and an intrinsic velocity of 1.1 ± 0.2 c. This large velocity 
potentially violates the speed of light and can be explained 
as either (i) the incorrect identification of the ejection date 
(i.e. not coincident with the X-ray flare), (ii) the jet had a 
much higher value for j3 cos 9 or (iii) an incorrect distance to 
the source has been estimated. As only one X-ray flare oc- 
curred during this period, it is unlikely that (i) can explain 
the high velocity. Also, as the same distance was used for 
all four knots (iii) does not explain the significantly higher 
velocity found in the first knot. Therefore, this effect is most 
likely due to (ii) - a change of the jet angle to the line of 
sight. 

Table [2] shows that the ejections in 2003 exhibited a 
much lower proper motion of app roxi mately 18 mas d -1 , as 
descri bed bv lFender et al.l l|l999h and lMirabel fc Rodriguez] 
l|l994h . The initial observation on 2003 March 06 re- 
vealed only a small extension associated with the start 
of an acti ve period. F o llowin g this, the VLBA obser- 
vation by iFuchs et al.l (|2003T l found evidence of only 
a highly compact jet on 2003 April 02 (MJD 52731), 



shortly before an X-ray flare, but MERLIN observations on 
2003 April 18 (MJD 52747), 13 days later, revealed the ejec- 
tion of a plasmon clearly associated with this event. The 
solitary X-ray flare on 2003 April 05 (MJD 52734) is there- 
fore the most probable progenitor of the observed plasmon 
and would make its proper motion 16.5±1 mas d _1 . Interest- 
ingly, GRS 1915±105 returned to the plateau state demon- 
strating that only a short (~ day) change in its X-ray state is 
required to eject material and the underlying emission from 
the plateau state had remained. This supports the idea that 
the knots were formed from internal shocks, rather than a 
cataclysmic destruction of the jet-disk model. It also demon- 
strated that the knots were formed soon after the X-ray 
spectral change. 

The period of activity in 2003 finally ended in a similar 
manner to the end of the 1997 event. MERLIN observations 
on 2003 June 15 (MJD 52805) showed another ejection of 
a possible superluminal plasmon. As this event ended with 
multiple X-ray flares an exact proper motion cannot be iden- 
tified from the MERLIN observations. 



4.4 Variations in polarisation 

GRS 1915±105 shows a large variation in the fraction of lin- 
ear polarisation in the core from less than 0.2 % to about 
14 %. During the ejection of material, the extended jet shows 
polarised components whilst the core appears depolarised. 
The electric field vectors of the jets are orientated along the 
jet axis if the rotation measure is low and the magnetic field 
is therefore perpendicular to the jet. This is suggestive of a 
tangled magnetic field compressed by a shock front travel- 
ling along the jet. 

The low levels of polarisation in the core can be ex- 
plained if the inner emitting jet passes through various 
depths of surrounding Faraday material. Differential degrees 
of Faraday rotation would cause the cancellation of linear 
polarisation, until the emitting material is far enough from 
the central region. 

As recently ejected components are likely to be in a 
much denser surrounding medium tha n those that were 
eject ed earlier and have moved away (|Miller- Jones et all 
l2007l ). one would expect the Faraday depth to fall as the 
components moves outward and depolarisation would be ex- 
pected to decrease and the frac tional polaris a tion would 
increase. This is in contrast to IFender et af] (|l999T ) who 
found the fraction of linear polarisation to quickly decay 
from ~ 13% to ~ 6% in the approaching jet. Addition- 
ally, iRodriguez et~ai"1 1 19951 ) reported no Faraday rotation 
in earlier VLA observations, showing no changing in po- 
sition angle between 5.0, 8.4 and 15 GHz, resulting in a 
RM< 50 rad m -2 . Observing no change in position angle 
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across the MERLIN sub-band is thus consistent with the 
VLA results and constrains RM< 600 rad m -2 (at a 3a 
limit). 

The multi-wavelength VLBA observations in 1997 were 
taken after a state change (from the plateau state) associ- 
ated with superluminal knots. Linearly polarised emissions 
were detected from the extended jets; however, no linear po- 
larisation was detected from the core. The results presented 
in Section [3] showed the linear polarisation was a higher 
fraction of the total jet emission at the higher frequency. 
Whilst the lack of linear polarisation from the core can be 
explained by intrinsic changes in the outflow (i.e. simply the 
absence of linearly polarised emission) , this cannot account 
for the changes with respect to frequency; changing levels 
of fractional polarisation with frequency is predicted by the 
presence of a Faraday rotating medium. These arguments 
therefore support the observed variation in polarisation as 
due to partial depolarisation by Faraday rotation. 



5 CONCLUSIONS 

GRS 1915+105 has been simultaneously observed in the ra- 
dio and X-rays in various different states over a ten year 
period. This comparative study between the X-ray and ra- 
dio variability shows, for the first time, both long term 
(~ weeks/months) and short lived (intra-day) variations. 
Daily monitoring shows periods of extended activity and 
X-ray spectral states that are related to the ejection of ma- 
terial. Radio activity over this period has shown to originate 
from the central 150 mas core, unless a short X-ray flare is 
observed that marks the ejection of a superluminal plasmon. 

The proper motion of ejecta has been calculated by 
identifying the time of ejection and exhibit no deceleration. 
Measured velocities were in agreement with previous obser- 
vations, and show a significant variation in velocity. Polari- 
sation measurements show strong linear polarisation varia- 
tion. Depolarisation in the core of the XRB is found during 
the ejection of a plasmon and the multi-wavelength observa- 
tions show an increase of linear polarisation with frequency, 
suggesting internal Faraday rotation effects. 

Observing GRS 1915+105 over a long period of time has 
shown the relativistic out-flow in different forms. If the ac- 
cretion conditions permit, collimated and steady jets persist 
over many weeks. Different apparent velocities of superlu- 
minal knots suggest either a precession of the jet angle or 
variation in their formation. Finally, the nature of the short- 
lived radio flares remains unknown; for example, it is difficult 
to explain why these bright, but short-lived events do not 
produce extended structure, in contrast to a flare during a 
state transition. 
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